Directed transport in driven optical lattices by gauge generation 
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We examine the dynamics of ultracold atoms held in optical lattice potentials. By controlling 
the switching of a periodic driving potential we show how a phase-induced renormalization of the 
intersite tunneling can be used to produce directed motion and control wavepacket spreading. We 
further show how this generation of a synthetic gauge potential can be used to split and recombinc 
wavepackets, providing an attractive route to implementing quantum computing tasks. 
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Introduction - In recent years enormous experimental 
progress has been made in creating and trapping ultra- 
cold atom gases When placed in an optical lattice 
potential these gases provide extremely clean and con- 
trollable implementations of interacting lattice systems, 
since parameters such as the interparticle interaction, the 
lattice depth and spacing are all readily tunable. Dissi- 
pation and decoherence effects are typically extremely 
weak, allowing the quantum coherent behaviour of these 
systems to be directly observed. 

In contrast to electronic systems, however, trapped 
atoms are uncharged, and so electric or magnetic fields 
cannot easily be used to produce or regulate transport. 
Due to their excellent coherence properties, one means 
of controlling the dynamics of the atoms is via quantum 
interference effects. A notable example is termed "coher- 
ent destruction of tunneling" (CDT) in which a periodic 
driving of the lattice causes the amplitude of the inter- 
site hopping to be renormalized ;2]. This renormaliza- 
tion has been seen directly in the expansion of trapped 
Bose-Einstein condensates [H-Il] , and has been used very 
recently to produce the fascinating phenomenon of "su- 
per Bloch oscillations" @, @] , and to induce the quantum 
phase transition 0] between a superfiuid and an insula- 
tor. 

In this Letter we show that as well as controlling the 
amplitude of the hopping, a periodic driving field can 
also be used to produce a tunneling phase, equivalent 
to a U(l) gauge potential. This gauge potential arises 
from the combined effect of the phase of the driving field, 
and the careful control of the switching condition. Al- 
though here we only consider one-dimensional lattices, 
the technique can also be similarly applied to create hop- 
ping phases in higher-dimensional systems. In this case 
the phases can be interpreted as Aharonov-Bohm phases 
picked up by a particle hopping from site to site, cor- 
responding to a synthetic magnetic field threading the 
lattice [§]. Other schemes have been devised to produce 
such gauge potentials in cold atom systems, including 
lattice rotations, state-dependent optical potentials 
or phase- imprinting [lfj. Our procedure, however, has 
an appealing simplicity requiring only the periodic vi- 
bration of the lattice potential, which is easily produced 



in experiment. We show how the driving can be used 
to control both the spreading and position of an initial 
wavepacket, and in particular, how a directed current of 
non- dispersing wavepackets can be induced. We shall 
also demonstrate how wavepackets can be split, guided, 
and recombined in a controllable and robust manner, ac- 
cessible to current experiment. 

Model - A gas of weakly interacting ultracold bosonic 
atoms can be described well by the Gross-Pitaevskii 
equation (GPE). When a sufficiently deep optical lattice 
potential is applied, the wavefunction will localize mainly 
in the potential minima defining the lattice sites, making 
it convenient to use a discretised form of the GPE 

l lH = " ( J ^' +1 + Jt ^'- 1 ) +g '^' |2 ^ +jV ^- (1) 

Here tpj denotes the system's wavefunction on lattice 
site j, and J describes the the tunneling amplitude be- 
tween nearest neighbor sites. Interactions between the 
bosons are given by a mean-field interaction, set by the 
nonlinearity parameter g. The time-dependent driving 
potential is assumed to rise linearly across the lattice 
[S 0, @, and has a time dependence given by V(t) = 
A + Ksin(ujt + <fi), where A is a static tilt of the lattice 
potential, and to and K are the frequency and amplitude 
respectively of the oscillating component. 

As an initial state we take a Gaussian wavepacket, 
ipj = N exp [— j 2 /2a 2 , + iOj\ , where oo is the initial width 
of the wavepacket measured in units of the lattice spac- 
ing, and N normalizes the wavefunction to unity. This 
choice of initial state mimics the experimental situation 
[H EL H| > hi which the condensate is prepared in a har- 
monic trap, and so typically has a Gaussian profile when 
transferred to the optical lattice. Note also that we ex- 
plicitly include a site-dependent phase term Oj in the 
wavefunction. 

Analysis - We first consider the non-interacting case 
(g = 0). The Hamiltonian describing the system ([I]) is 
then T-periodic in time, where T = 2ir /uj, and the natu- 
ral framework to describe its time evolution is Floquet 
theory. This reveals that in the high-frequency limit 
(lo > J) the time-dependent driven system can be de- 
scribed by an effective static Hamiltonian, whose param- 
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eters can be systematically evaluated by using perturba- 
tion theory ll| on the Floquet states. While these states 
are explicitly time dependent, being T-periodic functions, 
in the high-frequency limit their time variation is rather 
weak. This is the origin, for example, of the well-known 
negligible time-dependence of CDT [2J, as compared to 
the large oscillations observed at low drivin g fre quencies 
when dynamical localization instead occurs [T3| . 

To first-order the tunneling amplitudes are modified as 
J — > J e g = J(exp —i Jq V{t')dt' ), where (• • •) indicates 

a time-average over the driving period T. We restrict our- 
selves to considering the case of resonant driving, when 
A = nuj, which yields the result 



J e g gives the result 



J cS /J = e -iK/<-><*»+ e in(<h*/i)j n . 



(2) 



where J n is the nth Bessel function of the first kind. 
For the case of n = 0, a similar expression was obtained 
in Refll3l for a driven Bose-Hubbard model. There a 
ramped driving potential was used to adiabatically trans- 
form the ground state of the Hamiltonian to a stroboscop- 
ically current-carrying Floquet state. Here, however, we 
use the tunneling phase in a very different way, to control 
the non- equilibrium dynamics of an expanding atomic 
wavepacket. This, in conjunction with the weak time- 
dependence of the Floquet states in the high-frequency 
regime, means that the results we present are not strobo- 
scope, and have a negligible dependence on the moment 
within each driving period at which the system is mea- 
sured. 

From Eq.© we can immediately note the importance 
of the phase of the driving, <f>. For a cosinusoidal driv- 
ing ((f) = ±7r/2), the most frequently considered case 
in the literature, this result simplifies to yield J e g/J = 
(—l) n Jn (K/uj) - the well-known Bessel function renor- 
malization of tunneling found in CDT. For sinusoidal 
driving (<j> = 0), however, the tunneling additionally ac- 
quires a phase J e g/ J = exp [—i (K/uj — rnr/2)] J n (K/uj). 

It is natural to ask whether this tunneling-phase has 
physical implications, since it would appear that <fi can 
simply be eliminated by a shift of the time coordinate. 
It is important to note, however, that we consider the 
driving potential V(t) to be switched on at a specific 
moment t = 0, in common with experimental imple- 
mentations 0, Q. This gives the time origin, and thus 
the driving-phase, an unambiguous definition, and con- 
sequently <fi can indeed be of experimental relevance, as 
noted in RefQ This differs from many theoretical anal- 
yses [IH , in which the steady-state properties of a driven 
system are considered, and the driving is implicitly as- 
sumed to have been turned on at t — > — oo. In such cases 
the phase of the driving is indeed unimportant. 

The expansion of an initially Gaussian condensate in 
a periodically-driven lattice was analyzed in Reflla for 
real values of J e g. Extending this analysis to complex 



a(t) = a Jl + ($t[J e s]t/o%y 



(3) 



We thus see that the spreading of the wavepacket is gov- 
erned by the real component of J e g . For n = 0, for exam- 
ple, SJ[J c ff] = Jo (K/ui) cos (K/uj cos<p), and so as well 
as freezing at the "standard" CDT condition (when the 
Bessel function vanishes) , expansion is also suppressed at 
an additional set of values where cos (K/uj cos</>) = 0. 

As well as the expansion of the condensate, another 
useful experimental measurement is its center of mass 
motion. In the absence of driving our system has the 
standard spectrum of a non-interacting lattice model, 
Ek = — 2Jcosfc. When the system is driven, we can 
replace the energies Ek with quasienergies, obtained as 
solutions of the Floquet equation, to obtain the new dis- 
persion relation Sk — — 2| J e g| cos(fc — fco), where J e ff = 
| J e ff | exp[zfco]. The effect of the tunneling phase is thus 
not to alter the quasienergy spectrum of the system, but 
to displace the wavepacket to another point in the first 
Brillouin zone. In analogy with the familiar semiclassi- 
cal expression we can now define a mean group velocity, 
v g = dek/dk, where the average is taken over one period 
of the driving, to obtain the final result 



(4) 



We thus arrive at the rather elegant result that the 
two quantities most accessible to experiment - the 
wavepacket expansion and its center of mass motion - 
are directly related to the real and imaginary parts re- 
spectively of J e g. 

Directed transport - To verify these results we numer- 
ically simulate the model ([1]) for a 200 site lattice with 
no static tilt (A = 0), and take the onsite phases 9j to 
be constant. In FigQJ, we show the condensate's expan- 
sion for a cosinusoidal driving for several values of K/uj. 
These curves consist of an initial quadratic dependence 
on t followed by a linear ballistic expansion at long times 
[l6j ]. and clearly show how varying K/uj controls the con- 
densate spreading. Eq|3]can be used to extract the value 
of |3?[J e ff]| from these expansion curves, which we plot 
in FigQ]}. The expected Bessel function dependence of 
J e g is clearly seen, with the condensate expansion being 
frozen at the zeros of Jq (for K/uj = 2 AO, 5.52 . . .). How- 
ever, the corresponding expansion for sinusoidal driving 
shows an additional set of zeros at K/ui = n/2, 3%/2 . . ., 
in exact agreement with Eqf2]for n — 0. At these values 
of driving the suppression of the expansion arises from 
a very different cause; the tunneling phase displaces the 
wavepacket to fco = n/2 in the first Brillouin zone where 
the quasienergy bands have an inflexion, causing the ef- 
fective mass to diverge and so quenching the spreading 
of the wavepacket. 



3 



(a)' 








1 1 






-(c) 





Vs. 

\ \ 
\ \ 


1 

(b) 




•••• 




/ \ / 
(d) 





K/m 



FIG. 1: Response of a Gaussian wavepacket (oo = 4) to 
a periodic driving potential, lj — 16J. (a) Wavepacket 
expansion under cosinusoidal driving; from top to bottom 
K/lj = 0,1,2,2.4. For K/lj = 2.4 CDT occurs, and the ex- 
pansion is suppressed, (b) Real component of J e g, extracted 
from the expansion curves for cosinusoidal (black circles) and 
sinusoidal (red squares) driving. The curves show the the- 
oretical prediction obtained from Eq[3] (c) Displacement of 
a wavepacket under sinusoidal driving in units of the lattice 
spacing. For K/lj = and 2.4 no displacement occurs; oth- 
erwise it increases linearly with time, (d) As in (b), for the 
wavepacket velocity, given in units of cIl/T where dz, is the 
lattice spacing. Vertical blue arrows mark the driving pa- 
rameters K/lj = 7r/2 (dispersionless directed transport) and 
K/lj = 2.404 (complete suppression of dynamics). 



In FigfTJ; we show the motion of the center of 
mass of the condensate under sinusoidal driving. The 
wavepacket, initially at rest, begins to move at a con- 
stant rate, depending on K/uj. In FigQjL we plot the ve- 
locity corresponding to this displacement, and find that 
it agrees excellently with the predicted mean group ve- 
locity v g (Eq|4j. Under cosinusoidal driving, however, 
the velocity of the wavepacket is zero, also as predicted. 

We note that to obtain Eqs l3l4l we have assumed that 
the driving field can be turned on instantaneously. In ex- 
periment, of course, this idealized behaviour is not pos- 
sible. To check if our results are robust to this effect, we 
have included a ramp function in V(t) to describe the ef- 
fect of turning the field on from zero during a short, but 
finite, time interval. We find that as long as the ramp- 
time is sufficiently short, < 0.02T, very similar results are 
obtained. Experiments typically use driving frequencies 
of the order of kHz, which would thus demand ramp- 
times of ~ 10/iS, which are achievable. 

Directed motion - We can thus see that <f> can be used 
to cause an initially stationary wavepacket to move in a 
given direction with a precisely defined velocity, without 
requiring the spatial symmetry of the lattice to be bro- 
ken. Two values of K/lj are of particular interest, and are 
marked in FigQJ>,d. For K/uj = 2.404 (the first zero of 
Jo) the expansion of the initial wavepacket is suppressed, 
and its induced velocity is zero for all values of the driv- 
ing phase. This amounts to a complete suppression of the 
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FIG. 2: Density plot of the motion of a wavepacket under 
sinusoidal driving. Initially K/lj = n/2 to induce dispersion- 
less transport; K/lj is then set to 2.404 to freeze the motion, 
and finally to —n/2 to reverse it. The vertical dotted lines 
indicate the times at which K/lj is changed. 



dynamics of the condensate. However, at K/lj = tt/2 a 
wavepacket that is sinusoidally driven will not expand, 
but will have a non-zero velocity - a directed current of 
non-dispersive wavepackets. 

In Fig(5]we show the motion of such a wavepacket. Ini- 
tially we set K/uj = 7r/2 to induce motion. The driving 
is then tuned to K/lj — 2.404 to bring the wavepacket to 
a halt, and then to K/uj — —tt/2 to move the wavepacket 
in the opposite direction. It is clear that the spreading 
of the wavepacket is negligible, and that this technique 
indeed gives excellent control over the system. It is inter- 
esting to note that a similar form of control was reported 
in Ref 1 1 71 for an amplitude modulated lattice, instead of 
the phase-modulated lattice we consider. An important 
difference between the two cases, however, is that phase 
modulation does not require the presence of a static lat- 
tice tilt, since the effects also occur for n = 0, whereas 
amplitude modulation is limited to the case of resonant 
driving (n > 0). In addition, amplitude modulation does 
not produce CDT, the intersite tunneling depending lin- 
early on the driving amplitude instead of the Bessel func- 
tion dependence given in EqjSJ 



Wavepacket splitting - We now consider the effect 
of the onsite phases 6 a. It is well-known [l(| that im- 
printing a wavepacket with a uniform phase gradient 
0j+i — 0j = 9 has the effect of inducing motion of the 
center of mass, similar to the motion we have observed 
by manipulating the driving-phase <fi. By simulating the 
system with different values of 9, we have confirmed 
that the two phases combine, so that the net motion of 
the wavepacket actually depends on the phase difference 
<f) — 9. The driving field can thus be used to separate com- 
ponents of a wavepacket which possess different phase 
gradients. Let us consider the case of a superposition of 
a wavepacket with uniform phase (6 = 0) and one with tt- 
phase (9 = n). If the components have equal weight, the 
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FIG. 3: Density plot of a two-component wavepacket under 
sinusoidal driving. Initially K/uj — n/2 and the wavepacket 
splits apart. Its motion is then halted and reversed as before 
to bring about a collision. The interaction, g = 0.5J, causes 
the final state to be asymmetric. 

superposition will have the form ipj = N exp[— j 2 /2ctq] 
for j odd, and ipj = for j even (with no loss of gen- 
erality we can interchange the roles of the odd and even 
sites). Such a state can be prepared, for example, by 
patterned loading of a single uniform-phase condensate. 

Under sinusoidal driving, for K/oj = n/2, the compo- 
nent with uniform phase will move, without distortion, 
at a negative velocity, while the 7r-phase component will 
move identically but with a positive velocity. The initial 
wavepacket will thus split apart, as shown in Figj3] Tun- 
ing K/cj to 2.404 will bring each component to a stop, 
and then setting K/uj = —Tt/2 will bring the wavepackets 
together. For zero interaction (g — 0) the wavepackets 
will simply pass through each other. For small values of 
g the splitting process occurs as before, but during the 
collision the interaction causes the wavepackets to distort 
and produces a slightly asymmetric final state as shown. 



Incoherent expansion - An intriguing result seen in 
Refja is that unlike previous driven lattice experiments 
[H, S] , the wavepacket deformed under resonant driving, 
developing pronounced edges during its expansion. As 
a possible explanation of this effect we now look at the 
expansion of a phase incoherent wavepacket, by averaging 
over many realizations of random onsite phases 9j . The 
result in FigH^, is strikingly similar to the experimental 
observation. 



The phase effects we have discussed give a sim- 
ple explanation of this behavior. A phase incoherent 
wavepacket can be expressed as a mixture of many 
wavepackets, each with a random, but constant, phase 
gradient 9. Under the periodic driving each component 
will both develop a certain velocity and will spread, ac- 
cording to Eqs l3l4l As the components which spread 



FIG. 4: (a) Expansion of a phase incoherent wavepacket, ob- 
tained by averaging over 200 random realizations of the onsite 
phases. Driving parameters: n = 1, K/ui — 0.8. The density 
profile is shown in timesteps of 24T, showing the formation 
of a double peak structure as the wavepacket spreads, (b) 
Symbols show results obtained from simulations for the real 
component of the effective tunneling, 5ft[J e ff] (black circles), 
which controls the rate of the wavepacket spreading, and ve- 
locity (red squares), given in units oid^lT; dashed lines show 
the analytical results. The minimum values of 5i[J e g] align 
with the maximum velocity, and vice-versa. 

least have the highest velocity, while those that move 
more slowly spread more rapidly (see FigHJa), the ini- 
tial state will segregate with the rapidly moving com- 
ponents at the edges of the wavepacket remaining taller 
and narrower than the slower-moving components near 
the center. The edges of the wavepacket will move 
at the maximum speed, which for n = 1 is given by 
^max = \2J\{K/u)\. For the driving parameters used, 
our model predicts w m ax = 847 d^/s, where dj, is the 
lattice spacing, which compares well with the experimen- 
tally measured value of 869 <2l/s. We thus suggest that 
the unusual expansion seen in Refi is a consequence of 
the phase incoherence of the initial state, possibly arising 
from phase randomization produced by Wannier-Stark 
localization during the preparation of the system. 

Conclusions - We have shown how the phase of a driv- 
ing potential can be used to control the dynamics of an 
atomic wavepacket, both by regulating its rate of expan- 
sion, and by inducing a steady drift of its center of mass. 
Combining these effects allows the directed transport of 
non-dispersive wavepackets. Periodic driving also acts 
as a "prism" for the separation of different phase con- 
tributions within a wavepacket. This allows wavepack- 
ets to be divided and recombined, and also provides an 
appealing explanation for the unusual condensate expan- 
sion observed in Ref|6j. While these results have been 
obtained within mean field theory, probing the behavior 
of systems in the strongly-correlated regime remains an 
interesting subject for future research, holding out the 
enticing prospect of using these effects to generate and 
distribute entanglement in coherent lattice systems. We 
also note that since these directed currents require coher- 
ence across many lattice sites and driving cycles, their 
eventual decay may provide information on decoherence 
mechanisms. 
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